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' in Ising models with long-range interactions is investigated. We find the spinodal is 
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■ associated with the zeros of the partition function in four-dimensional complex tem- 
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perature/magnetic field space. The zeros approach the real temperature/magnetic 
field plane as the range of interaction increases. 

I. INTRODUCTION 



X 

i-j \ Mean-field treatments of fluids and Ising models yield metastable and unstable regions, 

separated by a well-defined line known as the spinodal As the spinodal is approached, 
the system shows phenomena similar to that at mean-field critical points. In particular, the 
isothermal susceptibility x diverges as a power law as the spinodal value of the magnetic 
field hs is approached from the metastable state P|, 

x^{K-h)-^/\ (1) 

where h is the (dimensionless) magnetic field. A mean-field Ising system can be realized in 
a well defined way by assuming an infinite-range interaction between the spins ^j. 

If the interaction range is long but finite, the system is no longer mean-field, but can be 
described as near-mean-field, and the spinodal singularity is replaced by a pseudospinodal. 
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which still has physical effects. As with apparent critical points in finite size systems, the 
susceptibility for a finite interaction range can be fit to a power law over a limited range 
of the scaling field {hs — h). For example, the susceptibility in the metastable state of 
a long but finite-range interaction Ising mode^appears to diverge over several decades in 
{hs — h) as the pseudospinodal is approached [21] • However, the divergence is suppressed if 
the pseudospinodal is approached too closely, indicating that the spinodal singularity has 
been smeared out 0. It also is found that the properties of the pseudospinodal converge 
rapidly with increasing interaction range to those predicted for the spinodal in mean-field 
theory Q]. 

Because the spinodal is a line of critical points, we expect that the Ising spinodal has 
properties similar to those of the Ising critical point. In particular, we expect that the 

n 

spinodal is related to the zeros of the partition function. Lee and Yang [5] showed that the 
singularity in the free energy for T less than T^, the critical temperature, arises from the 
presence of a real positive zero of the partition function in the thermodynamic limit. In 
finite systems there is no real positive zero for T < Tc. 

The zeros of the partition function for real temperature T lie on the imaginary axis of the 
complex magnetic field plane j^. For finite systems there is a gap in the distribution of the 
zeros around h = 0. This gap shrinks to zero for T < T^. as N ^ oo. The Lee- Yang relation 
between the zeros of the partition function and critical points is valid for all interaction 
ranges at the Ising critical point. This ideas were extended to complex temperature by 
Fisher Q. 

In order to generalize these ideas to spinodals in Ising models, we consider an Ising model 
in a magnetic field h. We will consider the "infinite-ran^e" Ising model in which each spin 
interacts with all other spins , and the Domb-Dalton jj version of the Ising model in which 
each spin interacts with its neighbors within a given interaction range R with a constant 
interaction. These models can be described by mean-field theory in the limit N ^ oo and 
R ^ oo (for infinite), respectively. Our main result is that the pseudospinodal is related 
to the zeros of the partition function in four- dimensional complex temperature/magnetic 
field space. In addition, the zeros approach the real temperature/magnetic field plane as 
the system becomes more mean-field. 

The structure of the paper is as follows. In Sec. O we consider the infinite-range Ising 
model and show both analytically and numerically that the zeros of the partition function 
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approach the real /3 and h plane as increases. We find a similar result in Sec. II III by 
estimating the partition function in the metastable state by using the Metropolis algorithm 
and the single histogram method. In Sec. IIVI we consider the Domb-Dalton version of the 
Ising model and estimate the partition function in the metastable state in the same way. 



II. THE INFINITE-RANGE ISING MODEL: ANALYTICAL APPROACH 

We first consider an Ising model in which every spin interacts with every other spin. We 
will refer to this model as the infinite-range Ising model j^, although the interaction range 
becomes infinite only in the limit — oo. The Hamiltonian is 

N 

H = -Jn ^ (TiUj - h^ai. (2) 

We need to rescale the interaction so that the total interaction energy seen by a given spin 
remains the same as is increased jH]. We will take 

* - l^V (3) 

This choice of Jn yields the Ising mean-field critical temperature Tc = 4 when oo Q], 
where we have chosen units such that J/k = 1. 

The exact density of states is easily calculated for this model and is given by 

A^' 

9{M) = ^T- rr, (4) 

nl[i\ — njl 

where n is the number of up spins. We have M = 2n — N and the total energy E = 
Jn{N — M'^)/2, where M is the magnetization. (In general, the density of states depends on 
both E and M, but because ii^ is a unique function of M in the infinite-range Ising model, 
we need only write g{M).) Hence the partition function Z can be expressed analytically for 
arbitrary complex inverse temperatures (3 and magnetic fields h: 

Z{p,ph) = ^ ^(M)e-^^e'3'^^. (5) 

allM 

To understand the nature of Z in the metastable state, imagine a simulation of an Ising 
model in equilibrium with a heat bath at inverse temperature (3 in the magnetic field h = 
hf) > 0. Because > 0, the magnetization values will be positive. Then we let h —h^. If 
/io is not too large, the system will remain in the metastable state for a reasonable amount 
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of time and sample positive values of M. Hence, to determine the zeros of Z associated with 
the pseudospinodal, we need to restrict the sum in Eq. (0) to magnetization values, M, that 
are representative of the metastable state. The following examples will illustrate the need 
for this restriction and the procedure for determining the zeros of Z. 

The notion of using a restricted partition function sum to describe the metastable state 



has a long history. Penrose and Lebowitz 10|, review such restricted partition functions 



and their properties. A more physical approach can be found in the discussion of the 
non- interacting droplet model by Langer [ij]. In this model fluctuations are restricted to 
non-interacting compact droplets of the stable phase occurring in the metastable phase. The 
partition function sum is restricted to droplets less than the critical size. This approximation 
is reasonable for low temperatures close to the coexistence curve. Langer showed that this 
restriction gives the same metastable state free energy as the analytic continuation of the 
stable state free energy in the same model. However, there are additional properties of the 
analytic continuation that do not appear in the restricted sum which are related to the decay 
of the metastable state rather than the description of the metastable state itself. 

We first consider = 4 and retain only the terms in the partition function sum that 
correspond to the two positive values of M. From Eqs. Q and (0), we have 

(3h) = 9{M)e-^^e^''^' = e^'^e"^^'^ + 4e-^^\ (6) 

M=2,4 

The subscript r denotes that the sum over M is restricted. If we let x = e"^^'*, the equation 
= is equivalent to, 

e^Px + 4 = 0, (7) 



and has the solution 



Ph= - \n2-i^+4p. (8) 



In general, we have four unknowns (the real and imaginary parts of f3 and h); Eq. (jH} 
yields two conditions. In the following we will fix 

Rep = 9/16, (9) 

which is equivalent to a temperature of T = |Tc. For this value of T, the value of the 
spinodal magnetic field is known to be hg ^ 1.2704 Equation ((Tj) then gives a line of 
zeros in complex [(3,(3h) space. However, if we are interested only in the zero closest to the 
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real (3, (3h plane, we need a fourth condition. This condition is found by requiring that the 
quantity, 

= (Im/3)2 + (Im/3/^)^ (10) 

be a minimum, which is equivalent to requiring that the leading zero of Z^, the zero closest 
to the real (3 and [3h plane, be as close to this plane as possible. If we let ?/ = Im/3 and use 
Eq. (jHJ, we can rewrite D"^ as 

D' = y' + {~\ + Ayf. (11) 
Because we want to be a minimum, we require 

^ = 2i/ + 2(-- + %)(4) = 0. (12) 

The solution is 

y = Im/3 = ^ ^ 0.3696, (13) 

and D = 0.38097. Note that ReT = Rep/{Rep^ + Im/?^) = 1.2417. We finally use Eq. © 
to obtain the value of complex h. The result is summarized in the first row of Table HI 

The solutions for = 9, 16, 100, and 1000, keeping all the positive M contributions to 
the partition function, also are shown in Table HI We see that although D becomes smaller 
as is increased, |Re /i| overshoots the mean- field value of /i^ ~ 1-27 (for j3 = 9/16). Hence, 
retaining all the positive M terms in the partition function allows the system to explore 
more than the metastable state, and we need to further restrict the sum over values of M. 
Physically we want to exclude values of M that would drive the system to the stable phase. 

What are the appropriate values of the magnetization that will keep the system in the 
metastable state? One way to determine these values is to look at P{M), the probability 
that the system has magnetization M for a particular value of f3 and h: 

P{M) = g{M)e-^^e^^^. (14) 

Figure^shows P{M) for a system of = 400 spins for h = —1.0 and (3 = 9/16. The negative 
magnetization values have a relatively high probability (because h < 0) and correspond to 
the stable phase. The positive values of magnetization have a much lower probability, and 
the peak at M = 360 corresponds to the most probable value of M in the metastable state. 
In between the peak at M = —400 and the peak at M = 360, P{M) has a minimum at 
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^min = 192 for this value of h and an inflection point at Mj = 298. We will only include 
values of the magnetization in the partition function sum that are greater than Mj. 

The reason for this choice of the cutoff has to do with the nature of metastability. We 
expect that P{M) oc exp{—(3F{M)), where F{M) is the metastable state free energy. We 
want to exclude from the partition function values of M that correspond to states which are 
not characteristic of equilibrium. In the infinite range model the configurations with Mmin < 
M < Mj are unstable in that the initial evolution of a fluctuation does not monotonically 
decay to the metastable well. This behavior translates into an initial growth of fluctuations 
rather than the monotonic decay expected in equilibrium. 

This behavior is a consequence of the fact that the free energy, which is proportional to 
log(P(M)), is not convex for these values of M. Obviously, this behavior holds for a range 
of values of M < Mmin as well. However, we can exclude all configurations with M < Mmin 
because they are in the stable free energy well and do not occur in the metastable state. 
Our particular choice of the cutoff is well defined, but is arbitrary to some extent as long as 
M is greater than Mj. 

To determine M/, we calculate the second derivative of P{M) as given in Eq. (fT^ . We 
find the value of M that satisfies 

d'PjM) _ -N 
dM^ ~ iN-M)iN + M)^^ ^ ^ 

Clearly there will be two inflection points (see Fig. P). We choose the one closest to the 
metastable state maximum of P{M). We find that the value of Mj is independent of h, 
which is consistent with the idea that the free energy for this system can be written in the 
Landau-Ginzburg form where the magnetic field appears only in a term linear in M. 
We now write the restricted partition function Zr{l3^(3h) as 

N 

Zr{P,Ph)= Cmx^"I\ (16) 

M=Mi 

The coefficients, Ca/, extend over a wide range of values and are as large as Wp'^^ for the 
values of that we considered. For this reason we computed Cm to arbitrary precision 
so as not to lose accuracy. The zeros of Z,,, which is a polynomial in a;, were found using 



MPSolve 



i3- 



For a given value of Im /3, we solve for the zeros of Z.,. in Eq. ()16p and find the value of x 
that corresponds to the leading zero, the zero that minimizes D in Eq. (fTUI). We repeat this 
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step for a range of values of Im (3 and determine numerically the value of Im (3 that yields 
the minimum value of D. The typical dependence of D on ImjS is shown in Fig. |21 From 
Fig. 121 we see that for N = 400, D is a minimum for ImjS ^ 0.035. Once we know this value 
of Im/?, we solve for h from the relation h = logx/(— 2/?). (The value of x was determined 
from the solution of Zr = 0.) 

We repeat the above steps for a range of values of and obtain D, Im/3, Keh, and 
Im h. Our results are summarized in Table |nj Note that Im h, Im (3, and D decrease as 
increases and |Re/i| approaches = 1.27. A plot of the zeros of Z for the infinite-range 
Ising model in the Imx, Rex plane is shown in Fig. IHl The values of D listed in Table ITTl 
are plotted as a function of A^ in Fig. El Because this log-log plot indicates a power-law 
dependence, we write 

D oc A^-^ (17) 

A least squares fit gives a = 0.659 ± 0.003. The estimate of the error is only statistical. 

Our numerical result for the exponent a can be understood by a simple scaling argument. 
In order for a mean-field approach, including the idea of a spinodal, to be a reasonable 
approximation, the system must satisfy the Ginzburg criterion, that is, the Ginzburg pa- 
rameter G must be much greater than unity. For the infinite-range Ising model, the Ginzburg 
criterion can be written as jl6| : 

^2 « 1' (18) 
where ^ is the correlation length, is the order parameter, and d is the spatial dimension. 
We have x ~ A/i~^^^ and ~ Ah^^'^ where Ah = hg — h, to obtain 

G = i'^Ah^l'^ » 1. (19) 

For the infinite-range model, A^ ~ Hence, 

G = NAh^/\ (20) 

up to a constant of order unity. The Ginzburg parameter G is a measure of how mean-field 
the system is for finite A^; the larger the value of G, the more mean-field the system is. If 
we keep G constant as we approach the spinodal, we see that 

Ah oc N-'^^\ (21) 

We present an argument for why G should be held constant in Sec. 



8 



From Eq. !^21\i we see that Ah and the distance D approach zero with an exponent 
a = 2/3, in good agreement with our numerical result. Because Ah in Eq. (PT|) could be 
associated with Imh, [hg — |Re/i|), or D, we note that Imh in Table ITTl goes to zero with 
the same exponent as in Eq. ()17|l . We also find that {hg — |Re/i|) oc N^^-^^. 



III. THE INFINITE-RANGE ISING MODEL: MONTE CARLO APPROACH 

In general, the partition function is not known analytically. However, we can use a 
Monte Carlo (MC) method to determine the density of states from which we can determine 
an estimate of the partition function. Such an app roach has been used to find the density of 
states for the nearest-neighbor Ising model [itIIisIi. In this way the leading partition function 
zeros at the critical point have been computed and the critical exponent u and corrections to 
scaling have been found with high precision Our goal is not to obtain precise estimates 
of the critical exponents near the pseudospinodal, but to show that the same simulations that 
show an apparent divergence in the susceptibility also yield an estimate for the leading zero 
of the partition function which behaves as expected as the system becomes more mean-field. 

To this end we use the Metropolis algorithm to equilibrate the system at temperature 
T = 16/9 and applied magnetic field h = ho for about 100 MC steps per spin. (The system 
equilibriates as quickly as 10 MC steps per spin depending on the strength of the field.) Then 
we fiip the magnetic field and compute the histogram H{E, M) from which we determine 
the densify of states g{E, M) and the partition function for complex /3 and h. We save the 
values of M after h —>■ —ho and run until M changes sign or for 5000 MC steps per spin, 
whichever comes first. We then throw away the first 20% of the data to ensure that the 
system is in metastable equilibrium and the last 20% of the data to ensure that we do not 
retain values of M that are too close to the stable state. The remaining 60% of the run is 
used to obtain H{E,M). We also omit any run whose lifetime in the metastable state is 
less than 100 MC steps per spin. Our results for H{E, M) are not sensitive to the choice 
of the minimum lifetime nor the percentage of each run that we use to estimate H{E, M). 
(Because ii^ is a function of M for the infinite-range Ising model, we need only to compute 
H{M). However, we need to compute H{E,M) in Sec. IV.) We averaged H{E,M) over 
approximately 5000 runs for each value of ho for a total of approximately 1.5 x 10^ MC steps 
per spin for a given value of h and A^. Our results for the susceptibility x are given in Fig. 
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As mentioned in Sec. HI x shows an apparent divergence with an mean-field exponent of 1/2 
until the pseudospinodal is approached too closely. 

Given the histo gram H{E,M) at Po = 9/16 and h = —ho, we use the usual single 
histogram method [isf and express the partition function for arbitrary (complex) /? and h 
as: 

ZU(3, (3h) = H{E, M)e(^°-'5)^e-(*''o-^'^)^. (22) 

E,M 

Note that we do not have to determine the lower cutoff for M because the Monte Carlo 
simulation only samples values of M while the system is in a metastable state (noted by the 
subscript m in Eq. ((221) )• As ho is increased for fixed N, the distance D initially decreases, 
but then begins to increase as hg is approached too closely. That is, D shows a minimum as 
a function of Hq. For each of value of we choose the value of Hq for which D is a minimum. 
A comparison of the histogram that was determined analytically in Sec. |n]and estimated by 
the Metropolis Monte Carlo algorithm shows similar qualitative behavior, except that the 
latter is approximately a Gaussian and extends to lower values of M (see Table lln|l . but 
with a smaller amplitude, [xx note change xx] 

Table UTTl shows our results for h, Im/3, and D for a range of values of A^ at the values of 
ho that minimize the distance D. A log-log plot of D versus A^ for the data in Table UTTl is 
shown in Fig. ^ a least squares fit gives a = 0.60 ± 0.03, which is consistent with the result 
obtained using the exact density of states (with a cutoff). Note that r, the average lifetime 
of the metastable state, for each value of A^ is approximately a constant at the value of ho 
that was chosen to minimize D (see Fig. Ej). We will use this fact to choose the value of ho 
for the long-range Ising model in Sec. IIVI 



IV. LONG-RANGE TWO-DIMENSIONAL ISING MODEL 

As discussed in Sec. HI the susceptibility of long-range Ising models in the metastable state 
shows an apparent divergence as the applied magnetic field is increased. In the following 
we show that this effect of a pseudospinodal in long-range Ising models is reflected in the 
behavior of the zeros of the partition function as a function of the complex temperature 
and magnetic fleld. We will show that as the interaction range R increases, the leading zero 
moves closer to the real plane. 

Following Refs. Q andO, we consider an Ising model such that each spin interacts with its 
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neighbors within a given interaction range R with a constant interaction J = A/q, where q 
is the number of interaction neighbors. (The factor 4 is included so that J = 1 for the usual 
Ising model on the square lattice.) If the thermodynamic limit is taken first the system 
is mean-field in the R ^ oo limit, and the system is described by Curie- Weiss theory jlll |. 
In this limit the metastable state ends at a spinodal point. The spinodal is a critical point 
and the susceptibility x diverges as in Eq. ([T)). 

We consider the Ising system on a square lattice with linear dimension L = 240 and 
= 57600. The interaction range R is defined such that a given spin interacts with any 
spin that is within a circle of radius R. The number of neighbors of a given spin is shown in 
the second column of Table IIVI This system is large enough [L is 10 times larger than the 
maximum value of R) for finite size effects to be minimal, but the finite size of the system 
implies that the zeros of the partition function must be complex for any range R. 

As in Sec. IIIH we equilibrate the system at inverse temperature P = 9/16 and applied 
magnetic field ho for 100 MC steps per spin. Then we flip the fleld and compute the histogram 
H{E, M) from which we determine the density of states g{E, M) and the partition function 
for complex (3 and h. The most time consuming part of our procedure is determining the 
appropriate value of Hq. Although we could determine by minimizing D as in Sec. IIIH 
we instead determined the mean lifetime r of the metastable state as a function of Hq. For 
small ho (away from the spinodal), r is greater than the duration of our runs which are 
2 X 10^ Monte Carlo steps per spin. As ho approaches hg, the lifetime begins to decrease. As 
mentioned in Sec. IIIH we found that the values of ho that minimize D are near hg and also 
yield approximately constant lifetimes for different values of A^. (This behavior was found 
for both the inflnite- and long-range Ising models.) We choose the largest value of for 
which the mean lifetime of the metastable state just begins to decrease below 2 x 10^ Monte 
Carlo steps per spin. This criterion for choosing ho is not as sensitive as minimizing D, but 
is much quicker although some additional error is introduced when ho is chosen this way. 
For a given ho, 10^ runs of 2 x 10^ MC steps per spin were done to determine the histogram 
H{E,M). 

The resulting values of the distance D as a function of the interaction range R are listed 
in Table IIVI and are plotted in Fig. |H| We flnd 

D oc (23) 
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The largest contribution to the above error estimate arises from the uncertainty in the values 
of the applied fields ho. Similarly, we find that the difference {hg — ho) varies with R as (see 
Fig.ini) 

hs-hoOC i?-l-97±0.06_ (24) 

The scaling behavior of {hg — |Re h\) is similar. The systematic error due to the uncertainty 
in ho is the largest contribution to the error estimates. 

The scaling behavior of D and hg — ho can be understood by a scaling argument similar 
to that given in Sec. |n] For a finite range system, the Ginzburg parameter is 

G = i?^A/i3/2-'^/4. (25) 
For two dimensions, Eq. ()25|) becomes 

G = R'^Ah. (26) 
Hence, if we keep the Ginzburg parameter fixed, we conclude that 

Ah oc R-\ (27) 

which is consistent with Eqs. and 



V. DISCUSSION 



We have shown that the pseudospinodal in Ising models has a well defined thermody- 
namic interpretation and can be associated with the leading zero of the partition function in 
complex temperature/magnetic field space in analogy with the behavior of the Ising critical 
point in finite systems. Our results for the nature of the approach of the leading zero of 
the partition function to the real temperature and magnetic field plane are consistent with 
simple scaling arguments. 

An essential ingredient in the scaling arguments was the condition that the Ginzburg 
parameter G was held constant as the system approached the pseudospinodal. As was seen 
in Sec. Ill, choosing the value of ho that minimizes the distance of the leading zero to the 
real temperature/magnetic field plane also leads to a metastable state lifetime that is found 
numerically to be constant. From nucleation theory near the spinodal, we know that the 
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lifetime of the metastable state, r, is given by 



r oc (28) 

where G is defined in Eq. For large G and R as well as small A/i (close to the 

pseudospinodal), it is easily shown that constant r implies constant G. To see this relation 
we simply replace i? by i? + bR and A/i by A/i + 6h in Eq. (|28|) and demand that r remain 
constant. For 6R and small, constant r implies constant G. Because G is constant, the 
scaling arguments of Sees. Ill and IV follow. 

The relation between the zeroes of the partition function and the spinodal provides a 
mathematical foundation for the notion of a pseudospinodal and clarifies the extent to which 
spinodals act like critical points. It also provides a possible route by which pseudospinodals 
in supercooled liquids can be characterized 
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Tables 



TABLE I: Results for the infinite-range Ising model if all positive values of M are retained in Z. 
For larger A'', |Re/i| overshoots hg = 1.27 and goes to zero as N is increased still further. 



N Im/3 |Re/i| Im/i D 

4 0.3696 1.8577 1.3849 0.3810 

9 0.1704 1.1987 0.4782 0.1823 

16 0.1070 1.1359 0.2926 0.1153 

100 0.0145 0.8854 0.0366 0.0164 

1000 0.0014 0.8335 0.0035 0.0016 



TABLE II: Values of Mj/N, \Reh\, Im h, Im P and the distance D for increasing values of A'" for the 
infinite-range Ising model. As explained in the text, the inflection point of P(M) determines Mj, 
the cutoff for M. Note that |Re/t| approaches hg = 1.27 and Mj/N approaches mg = 0.745356. 
For AT = 4000000, Mi/N = 0.745356. 



N Mi/N Imp \Reh\ Imh D 

100 0.7400 0.086 1.1601 0.2257 0.0902 

400 0.7450 0.035 1.2125 0.0949 0.0367 

800 0.7450 0.022 1.2322 0.0603 0.0230 

1200 0.7450 0.017 1.2407 0.0453 0.0176 

1600 0.7450 0.014 1.2456 0.0375 0.0145 

2400 0.7458 0.011 1.2508 0.0280 0.0112 
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TABLE III: Results from Monte Carlo simulations of the infinite-range Ising model. The simula- 
tions were done in the applied field —h^ and at the inverse temperature 9/16. The values of ho for 
each value of N were chosen so that the distance D to the real P-Ph plane is a minimum. As noted 
in the text, this criterion for ho also yields metastable state lifetimes that are roughly constant for 
the different values of N. The values of Mcut represent the smallest values of M that were sampled 
in the metastable state. 



N 


ho 




r 


Im/3 


\Reh\ 


Imh 


D 


128 


0.9 


0.23 


4808 


0.0144 


1.181 


0.0165 


0.0163 


400 


1.0 


0.44 


> 5000 


0.0070 


1.196 


0.0095 


0.0076 


800 


1.1 


0.56 


> 5000 


0.0050 


1.229 


0.0063 


0.0056 


2400 


1.205 


0.59 


4995 


0.0026 


1.256 


0.0043 


0.0027 


4000 


1.226 


0.63 


4992 


0.0021 


1.264 


0.0033 


0.0022 


8000 


1.246 


0.65 


4908 


0.0016 


1.269 


0.0024 


0.0017 



TABLE IV: Summary of results for the Ising model with interaction range R on the square lattice 
with linear dimension L = 240. The number of neighbors q within an interaction range R is given 
in the second column. The value of is determined for each R by choosing the lifetime of the 
metastable state to be approximately 2 x 10^ MC steps per spin. The duration of each run was 

2 X 10^ MC steps per spin and each run was repeated 10"^ times. Because we did not use the first 
and last 20% of each run, the total number of MC steps per spin for each value of R was 1.2 x 10^. 



R 


q 


ho 


Im/3 


|Re^| 


Im h 


D 


6 


112 


0.95 


0.0036 


0.992 


0.0084 


0.0038 


8 


196 


1.05 


0.0035 


1.011 


0.0071 


0.0035 


12 


440 


1.18 


0.003 


1.199 


0.0078 


0.0031 


15 


708 


1.215 


0.0025 


1.230 


0.0055 


0.0026 


18 


1008 


1.235 


0.0013 


1.252 


0.0039 


0.0014 


24 


1792 


1.248 


0.00095 


1.259 


0.0022 


0.0010 



Figures 
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FIG. 1: The probability of the magnetization, P{M), obtained from Eq. 1)14(1 for the infinite-range 
Ising model with N = 400, h = —1.0, and (3 = 9/16. The values of P{M) are plotted on a 
logj^Q-linear scale because of the dominance of the negative values of M. We include terms in the 
partition function sum over M from the inflection point of P{M) at Mj = 298 to M = 400 in the 
partition function. The inset shows the region where the inflection point is (see arrow). 
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FIG. 2: The value of D, a measure of the distance of the leading zero to the real [3 and [3h plane, 
versus Im/3 for the infinite-range Ising model for = 400. The minimum distance to the real axes 
occurs at Im [3 ~ 0.035 for this value of N . 
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FIG. 3: The imaginary versus the real part of the zeros of the partition function plotted in terms 
of the variable x = e~^^'* for = 400 (empty circles) and 1600 (filled circles) for the infinite-range 
Ising model. The zeros were obtained by the analytical method described in Sec. ^ using Im/3 
listed in Table ITll 
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FIG. 4: Log-log (base 10) plot of D, the distance of the leading zero of the partition function to 
the real axes (see Eq. (fTUIl l. versus the system size N, using the analytical approach described 
in Sec. ^ The slope is —0.659 it 0.003. The data is from Table ITll 
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FIG. 5: Log-log (base 10) plot of the isothermal susceptibility x as a function of {hg — ho) for 
the infinite-range Ising model with = 10000. The system was equilibrated using the Metropolis 
algorithm at a temperature T = 16/9 and applied field h = Hq. Then the field was flipped and the 
values of the magnetization were sampled in the metastable state. Note that x exhibits mean-field 
behavior over about 2 decades and the apparent divergence of x at the spinodal field hg = 1.27 is 
rounded off when (hg — ho) becomes too small. This behavior is an example of the influence of a 
pseudospinodal. The straight line with a slope of 1/2 is drawn as a guide to the eye. 
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FIG. 6: Log-log (base 10) plot of D, the distance of the leading zero of the partition function 
to the real axes, versus the system size N for the infinite-range Ising model obtained from Monte 
Carlo simulations. A least squares fit gives a slope of —0.60 it 0.03. The data is from Table ITTII 
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FIG. 7: The lifetime r of the metastable state in the infinite-range Ising model as a function of 
the applied field /iq for N = 128. The behavior of r for the long-range Ising model considered in 
Sec. IIVI is similar, and for the latter we choose /iq to be the field at which r(/io) just begins to 
decrease sharply. 
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FIG. 8: Log-log (base 10) plot of the distance D to the real axes versus the interaction range 
R for the two-dimensional Ising model. Note that the leading zero of the partition function moves 
closer to the real axes as the system becomes more mean-field. A least squares fit gives a slope of 
-1.8 ± 0.3. The data is from Table ITVl 
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FIG. 9: Log-log (base 10) plot of the difference hs — ho as a function of the interaction range R. 
A least squares fit gives a slope of —2.01 ± 0.07. The data is from Table HVl 



